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1 6.  Abstract 

\This  study  demonstrates  the  utility  of  semi-empirical  Loran-C 
time1  difference  (TD)  grid  calibration  techniques.  Theory  is  employed 
to  determine  the  functional  dependence  of  TDs  on  range  and  bearing 
from  the  Loran-C  chain  stations.  TD  measurement  data  are  then  util¬ 
ized  to  calibrate  the  unknown  coefficients  incorporated  in  the  semi- 
empirical  TD  model.  A  semi-empirical  model  is  presented  for  the 
West  Coast  Loran-C  chain  where  at-sea  TD  measurement  data  in  Southern 
California  revealed  large  discrepancies  between  U.S.  Coast  Guard 
original  predictions  and  measurements.  A  significant  reduction  in 
the  TD  errors  is  achieved  with  the  semi-empirically-calibrated  model 
relative  to  the  U.S.  Coast  Guard  original  grid.  The  accuracy  of  the 
calibrated  West  Coast  Loran-C  grid  is  further  evaluated  by  comparing 
the  calibrated  grid  with  measurements  which  are  not  used  in  model 
calibration.  Results  are  also  presented  which  show  the  sensitivity 
of  the  model  accuracy  to  the  quantity  and  distribution  of  measurement 
data  used  to  calibrate  the  model.  Guidelines  are  formulated  to  aid 
in  the  design  of  data  collection  requirements  for  future  semi-empiri¬ 
cal  grid  calibration  efforts.!' 
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1. 


INTRODUCTION 


1 . 1  BACKGROUND 

Semi-empirical  Loran-C  time  difference  (TD)  grid  cali¬ 
bration  techniques  have  been  successfully  employed  to  develop 
an  accurate  (approximately  100  nsec  rms)  calibrated  grid  for 
the  St.  Marys  River  Loran-C  chain  where  groundwave  signal  paths 
exhibited  "nearly  homogeneous"  signal  propagation  properties 
(Ref.  1).  The  study  reported  herein  extends  the  utility  of 
semi-empirical  techniques  to  the  development  of  a  calibrated 
Loran-C  grid  for  the  Coastal  Confluence  Zone  (CCZ)  where  sig¬ 
nal  paths  must  be  considered  as  mixed,  i.e.,  part  land  and 
part  sea  water.  In  particular,  semi-empirical  TD  grid  cali¬ 
bration  techniques  are  applied  to  the  West  Coast  Loran-C  chain 
where  the  U.S.  Coast  Guard  predictions  are  reported  to  result 
in  large  charting  errors,  especially  in  the  CCZ  between  Los 
Angeles  and  San  Diego. 


1.2  WEST  COAST  LORAN-C  CHAIN  CONFIGURATION 

The  U.S.  West  Coast  Loran-C  chain  is  a  long  base¬ 
line  chain  consisting  of  four  transmitters  (Table  1.2-1)  and 

•k 

two  monitor  Sites  as  illustrated  in  Fig.  1.2-1.  Time- 
sequenced  groups  of  pulsed  100  kHz  radiowave  signals  are 
transmitted  by  the  four  stations.  The  difference  in  time- 
of-arrival  (TOA)  of  signals  from  two  of  the  stations  is  a 


*The  time  difference  measurement  data  provided  for  this  study 
were  collected  with  the  monitor  at  Point  Pinos  controlling 
the  chain. 
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TABLE  1.2-1 

WEST  COAST  LORAN-C  CHAIN,  RATE  9940(SS6),  CHARACTERISTICS 


STATION 

COORDINATES* 

STATION 

FUNCTION 

EMISSION 
DELAY  & 
BASELINE 
LENGTH 
(Msec) 

RADIATED 

PEAK 

POWER 

(kW) 

LATITUDE  & 
LONGITUDE 

Fallon, . 
Nevada 

39°  33'  06 . 62"N 
118°  49’  56. 37"W 

M 

Master 

400 

George , 
Washington 

47°  03'  47 . 99"N 
119°  44’  39.5"  W 

W 

Secondary 

13796.90 

2796.90 

1200 

Middletown, 

California 

38°  46'  56.99"N 
122°  29’  44. 53"W 

X 

Secondary 

28094.50 

1094.50 

400 

Searchlight , 
Nevada 

35°  19’  18.18N 
114°  48’  17.43’’W 

Y 

Secondary 

41967.30 

1967.30 

500 

*Based  on  World  Geodetic  System  (WGS)  *  1972  Datum. 


measure  of  the  difference  in  distance  from  the  point  of  ob¬ 
servation  to  each  of  the  two  stations.  The  locus  of  all  points 
having  the  same  observed  difference  in  distance  to  a  pair  of 
stations  is  a  hyperbola,  referred  to  as  a  line-of-position 
(LOP).  The  intersection  of  two  or  more  LOPs  defines  a  user's 
position  when  compared  to  a  chart  containing  a  grid  of  cali- 
\  brated  LOPs. 


,3  os™ 


The  objectives  of  the  study  reported  herein  are  to: 


Develop  a  semi-empirical  calibration 
technique  (or  model)  for  the  West  Coast 
Loran-C  chain  TD  grid  with  an  accuracy 
equal  to  or  better  than  the  current  U.S. 
Coast  Guard  grid  calibration  procedures 
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•  Quantify  the  semi-empirical  calibration 
technique  accuracy  sensitivity  in  terms 
of  calibration  data  requirements. 

Based  on  the  sensitivity  analysis  results,  data  collection 
guidelines  are  formulated  for  future  semi-empirical  grid  cali¬ 
bration  efforts.  The  outputs  of  the  study  include  a  TD  grid 
calibration  model  and  computed  TDs  at  a  number  of  at-sea  data 
points  which  are  not  used  to  calibrate  the  grid. 

1.4  TECHNICAL  APPROACH 

Loran-C  TD  measurement  data  are  analyzed  to  identify 
dominant  propagation  parameter-dependent  trends  in  the  data. 
These  trends  and  theory  are  used  to  establish  separate  poly¬ 
nomial  structures  for  the  land  and  sea  water  signal  propaga¬ 
tion  phase  delay  models.  Phase  delay  over  a  mixed  path  is 
then  obtained  by  applying  Millington's  method  (Ref.  5)  which 
empirically  combines  land  and  sea  phase  delays.  Data  are  then 
used  in  a  Kalman  estimation  algorithm  (Ref.  3)  to  compute  the 
uncertain  coefficients  of  the  land  and  sea  models  of  the  TD 
grid  calibration  algorithm.  The  calibrated  algorithm  is  used 
to  compute  TDs  at  each  data  site  and  the  TD  residuals  (differ¬ 
ence  between  measurement  and  calibrated  TDs)  are  examined. 
Adjustments  are  then  made  to  the  TD  model  structure  in  an 
attempt  to  further  reduce  the  residuals.  This  process  of  ad¬ 
justing  the  model  structure  is  repeated  until  the  residuals 
agree  with  the  expected  theoretical  covariance  associated  with 
the  TD  model.  The  model  which  exhibits  the  "best"  performance 
is  then  selected  from  the  several  candidate  models  as  the  West 
Coast  TD  grid  model.  The  performance  of  the  West  Coast  model 
is  further  evaluated  against  measurements  not  used  in  calibrat¬ 
ing  the  model. 


Next,  the  sensitivity  of  the  selected  West  Coast  TD 
grid  model  is  assessed  as  a  function  of  density  and  quantity 
of  the  calibration  data.  This  is  accomplished  by  recalibrat¬ 
ing  the  coefficients  of  the  model  with  a  subset  of  the  data. 

The  resulting  TD  residuals  are  evaluated  for  a  number  of  uni¬ 
formly  distributed  data  densities  as  well  as  several  clustered 
data  sets.  Sensitivity  analysis  results  are  used  to  develop 
general  data  collection  requirements  in  terms  of  expected  accu¬ 
racy  of  the  semi-empirically  calibrated  TD  grid. 


1 . 5  REPORT  OVERVIEW 

Theoretical  basis  of  the  semi-empirical  TD  grid  cali¬ 
bration  technique  is  described  in  Chapter  2.  The  measurement 
data  analysis  results  used  to  identify  significant  trends  in 
the  data  are  presented  in  Chapter  3.  The  grid  calibration 
procedure,  semi-empirical  models  and  the  resulting  calibrated 
model  algorithm  for  the  West  Coast  chain  along  with  the  model 
performance  are  presented  in  Chapter  4.  The  sensitivity  of 
the  calibrated  model  accuracy  to  the  quantity  and  distribution 
of  data  is  assessed  in  Chapter  5.  Conclusions  and  recommenda¬ 
tions  are  presented  in  Chapter  6.  Calibrated  grid  model  algo¬ 
rithm  details  and  computed  TDs  at  the  calibration  data  sites 
appear  in  Appendix  A. 


2. 


THEORETICAL  BASIS 


2 . 1  INTRODUCTION 

Loran-C  groundwave  navigation  position  fix  accuracy 
is  primarily  dependent  on  a  chart  maker's  ability  to  accurately 
compute  (or  determine)  the  expected  TOA,  or  TD  which  is  the 
difference  between  two  TOAs,  of  received  groundwave  signals 
from  Loran-C  transmitting  stations  for  comparison  with  Loran-C 
receiving  equipment  outputs.  Mathematical  algorithms  are  used 
in  automatic  equipment  and  manual  procedures  depend  on  the 
availability  of  accurate  charts  of  Loran-C  LOPs .  In  either 
case,  accurate  knowledge  of  the  signal  phase  delay  which  es¬ 
tablishes  the  true  LOP  is  critical  to  Loran-C  position  fix 
accuracy. 


The  phase  delay  of  a  groundwave  signal  is  generally 
expressed  as 


4>  =  T  +  SF 

=  |  R  +  SF  (2.1-1) 

where  n  is  the  surface  refractive  index,  c  is  the  speed  of 
light  in  a  vacuum,  R  is  the  range  to  the  transmitting  station, 
and  SF  is  the  phase  of  the  secondary  factor  (Ref.  4).  The 
primary  signal  phase  delay,  T,  (also  referred  to  as  the  pri¬ 
mary  phase  delay,  or  primary  travel  time  in  Ref.  2)  is  the 


*Phase  delay,  phase,  propagation  delay,  travel  time  and  time 
delay  are  used  interchangeably  throughout  this  report,  and 
are  expressed  in  units  of  time. 


computed  travel  time  of  the  Loran-C  pulse  over  a  distance  equal 
to  the  transmitter- to-receiver  great  circle  path  length,  ac¬ 
counting  only  for  the  velocity  of  light  and  the  index  of  re¬ 
fraction  of  the  atmosphere.  The  phase  of  the  secondary  factor 
(referred  to  as  the  secondary  phase  delay  in  this  report)  is 
characterized  in  Ref.  4  as  a  correction  to  the  primary  phase 
delay  to  account  for  the  phase  delay  due  to  signal  propagation 
over  the  inhomogeneous  and  irregular  surface  of  the  earth.  In 
the  groundwave  phase  delay  equation,  Eq.  2.1-1,  T  is  the  domi¬ 
nant  term  and  involves  well-known  parameters.  The  SF  is  usual¬ 
ly  an  order  of  magnitude  smaller  than  T  but  significantly  more 
complex  to  determine  due  to  the  inhomogeneous  electrical  prop¬ 
erties  and  irregularities  of  the  earth's  surface. 

2.2  SF  COMPUTATION  TECHNIQUES 

A  number  of  analytical  and  empirical  SF  computation 
techniques  (or  models)  have  been  reported  in  the  literature 
(Refs.  4  through  7).  The  most  commonly-used  models  include: 

•  Homogeneous/smooth  path  model 

•  Mixed  Path  --  Millington's  empirical 
method 

•  Inhomogeneous  Path  --  Integral  Equation 
approach . 

*The  secondary  phase  delay  is  equal  to  the  sum  of  the  Secondary 
Phase  Factor  and  the  Additional  Secondary  Factor,  as  defined 
in  Ref.  2.  The  Secondary  Phase  Factor  is  a  correction  to  the 
primary  phase  delay  on  the  presumption  the  path  is  entirely 
sea  water.  The  Additional  Secondary  Factor  is  a  correction 
to  the  Secondary  Phase  Factor  which  accounts  for  an  inhomo¬ 
geneous  earth's  surface. 
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The  homogeneous/smooth  path  model  (Ref.  4)  is  useful  for  SF 
computations  over  a  homogeneous  (i.e.,  uniform  electrical  prop¬ 
erties)  signal  propagation  path  along  a  smooth  earth,  such  as 
an  all-sea  water  path. 

Millington's  empirical  approach  (Ref.  5)  is  useful 
for  computing  the  SF  over  a  mixed  (multiple-homogeneous  seg¬ 
ment)  path.  This  approach  empirically  combines  SFs  of  various 
homogeneous  segments  of  a  mixed  path.  For  example,  for  a  two 
segment  (land  and  sea  water)  mixed  path,  as  shown  in  Fig.  2.2-1, 
Millington's  formula  for  the  SF  over  the  mixed  path  of  length 

TL+TS  1S 

SF(aL,  cts ,  Tl  +  Tg)  =  j  [*A  +  <t»B]  (2.2-1) 

where 

♦A  =  SF(aL>TL)  +  [SF(as,TL+Ts)  -  SF(<jg,TL)]  (2.2-2) 

4>g  =  SF(0g,Tg)  +  [SF(oL,TL+Tg)  -  SF(aL,Tg)]  (2.2-3) 

cr.  is  the  conductivity  of  the  homogeneous  land 
segment  of  the  path 

Og  is  the  conductivity  of  the  sea  water 
segment  of  the  path 

SF(a,T)  is  the  SF  for  a  homogeneous  path  of 
conductivity  o  and  length  T 
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Figure  2.2-1  Two-Segment  Mixed  Path 
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The  accuracy  of  Millington’s  approach  is  known  to  be  good  (Ref 
8)  provided  that  reasonably  accurate  estimates  of  the  homogen¬ 
eous  segment  SFs  are  available. 

Millington's  approach  has  also  been  used  by  the  U.S. 
Coast  Guard  to  generate  "effective"  conductivity  maps  for  a 
number  of  operational  Loran-C  chains.  The  map  generation  pro¬ 
cedure  is  to  use  baseline  and  selected  land  TD  measurement 
data  to  estimate  the  effective  conductivity  along  each  of  the 
homogeneous  segments  of  the  mixed  paths  included  in  the  data. 
Conductivities  along  measurement  data  paths  are  then  adjusted 
to  vary  until  predictions  and  data  are  in  agreement  within  100 
nsec.  Effective  conductivity  values  are  extended  to  regions 
not  traversed  by  data  paths  by  defining  geological  structures 
which  are  expected  to  have  similar  conductivities. 

When  the  propagation  path  is  inhomogeneous  and  the 
terrain  is  irregular,  such  that  it  cannot  be  modeled  satisfac¬ 
torily  by  either  the  homogeneous  path  or  Millington's  mixed 
path  model,  a  more  sophisticated  and  complicated  integral  equa 
tion  model  (Ref.  6)  can  be  used.  However,  the  numerical  solu¬ 
tion  of  the  integral  equation  is  generally  expensive  and  cum¬ 
bersome  except  for  simple  terrain  irregularities  and  requires 
a  relatively  large  computer  storage  capability  to  process  all 
of  the  physiographic  data  characterizing  the  path. 

In  summary,  analytical  prediction  models  are  useful 
if  the  modeled  propagation  path  scenario  closely  approximates 
the  "real-world"  scenario  and  if  the  propagation  path  parame¬ 
ters  are  known.  Usually,  the  real-world  signal  propagation 
medium  of  interest  is  far  too  inhomogeneous  and  irregular  to 
be  easily  idealized-  Additionally,  the  required  propagation 
path  parameter  values  are  rarely  known  with  the  required  pre¬ 
cision. 
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2.3  SEMI -EMPIRICAL  SF  MODEL 

The  approach  taken  herein  is  to  employ  semi -empirical 
grid  calibration  techniques,  similar  to  those  used  for  calibrat¬ 
ing  the  St.  Marys  River  Loran-C  chain  (Ref.  1).  The  "physics" 
of  the  propagation  medium  are  used  to  establish  a  functional 
form  of  the  signal  phase  delay  model  and  measurement  data  are 
used  to  calibrate  the  (uncertain)  coefficients  of  the  model. 


A  generalized  semi -empirical  polynomial  functional 
form  for  the  SF  of  the  station  is  given  by 

K2  L 

SF<W  =  kgK  Ak  Tj  +  [Cj*  sin  «j  +  Dj*  cos 

1  (2.2-4) 

j  =  secondary  (w,  x,  or  y)  or  master  (m) 
station 

n  th 

Tj  =  —  Rj  =  j  station- to-user  primary  phase  delay 

Rj  =  j**1  station- to-user  great-circle  path  length 

Pj  =  user  path  bearing  angle  at  the  station 

,  K2  and  L  are  positive  integers 

C..  and  D..  are  the  station-dependent  coefficients 
J  J  of  harmonic  terms  in  the  model 

A^  is  the  range-dependent  coefficient  of  the 

model  which  may  in  general  be  station-dependent. 

The  semi-empirical  model  can  be  made  as  complex  as  desired  and 
will  approach  the  theoretical  model  in  the  limit.  However, 
increased  complexity  requires  estimating  an  increased  number 
of  uncertain  coefficients  in  the  model,  which  in  turn  increases 
the  amount  of  measurement  data  required.  Since  the  primary 


sfj  * 


where 
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purpose  for  developing  a  grid  calibration  model  is  to  reduce 
the  amount  of  measurement  data  required  to  establish  a  Loran-C 
grid,  a  compromise  must  be  made  between  model  complexity  and 
measurement  data  requirements.  For  example,  for  Loran-C  signal 
propagation  over  a  homogeneous/smooth  propagation  region  such 
as  an  all-sea  water  path,  the  SF  behavior  is  expected  to  be 
isotropic,  i.e.,  independent  of  both  path  bearing  angle  and 
station  location.  Therefore,  for  this  case,  the  semi-empirical 
SF  model  for  the  j*”  station  would  be  of  the  form 

k2 

SF.  =  SF(T . )  =  £  A.  (2.2-5) 

J  J  k=-K1  K  J 

where  coefficient  is  station-independent,  i.e.,  it  has  the 
same  numerical  value  for  every  station  of  the  chain. 


2.4  TD  GRID  CALIBRATION  EQUATION 


The  true  time  difference  (TD.)  between  the  time-of- 
arrival  (TOA.)  of  a  groundwave  signal  from  the  i  (=  w,  x  or 
y)  secondary  station  and  the  time-of-arrival  (TOAffl)  from  the 
master  station  (m),  is 


Time  Difference 


Time-of -Arrival 
From  Secondary 


TOA. 

1 


The  time-of-arrivals  can  be  expressed  as 


Time-of -Arrival 
From  Master 


TOA_ 

m 


(2.3-1) 


TOAi  =  Tt  +  SFt  +  EDi 


TOA_  =  T  +  SF 
m  m  m 


(2.3-2) 

(2.3-3) 
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where  ED.  is  the  true  emission  delay;  for  the  West  Coast  chain 
i  th 

it  is  equal  to  the  coding  delay  of  the  i  secondary  station 

of  the  chain  plus  the  true  baseline  length.  (Note,  the  published 

values  of  ED^  for  the  West  Coast  chain  are  as  given  in  Table 

A.l-1.).  Combining  Eqs .  2.3-1  through  2.3-3,  the  true  TD  is 

given  by 


TDi  =  (Ti"Tm)  +  (SFi_SFm)  +  EDi 


(2.3-4) 


The  semi-empirical  grid  calibration  model  developed 
herein  uses  land  and  sea  TD  measurement  data  to  calibrate  the 
model.  These  measurements  are  corrupted  by  measurement  noise 
including  position  reference  errors,  and  are  related  to  the 
true  TD  by 


Measured 

tda 


True  +  Measurement 
TD^  Noise 


* 

2i 


=  TD± 


vi 


(2.3-5) 


Upon  substituting  Eq.  2.3-4  into  Eq .  2.3-5,  the  measured  TD  is 


zi  *  <W  +  <SFi-SFm>  +  EDi  +  vi 


(2.3-6) 


In  subsequent  discussions  of  the  TD  data  quality  and 
model  calibration  procedure,  it  is  convenient  to  transform  the 
measured  TD^  into  an  "Adjusted  TD/’  which  is  defined  as 


ATDj  =  =  2*  -  <Tt  -  T0)  -  ED1  (2.3-7) 

where  ED^  is  the  published  constant  emission  delay  (see  Table 
1.1-1)  implemented  at  the  secondary  station.  Substituting  Eq. 
2.3-6  into  Eq.  2.3-7  gives 


(2.3-8) 


*  <SFt  -  SFn>  *  AED*  ♦  vj 

*  <SFi  -  SFm>  *  vi 

» 

where  (=  AED^  +  v^)  is  the  total  measurement  error,  and 

AED.  is  the  difference  between  the  true  and  published  emission 
1  th 

delay  for  the  i  secondary  station.  Equation  2.3-8  is  used 
in  Chapter  4  for  calibrating  the  SF  model.  Equations  2.2-4 
and  2.2-5  are  used  to  provide  the  basic  functional  structure 
for  the  SF  associated  with  all-land  and  all-sea  water  paths, 
respectively.  For  a  mixed  path,  the  SF  is  computed  by  using 
Millington's  empirical  method,  Eqs.  2.2-1  through  2.2-3.  There¬ 
fore,  the  task  of  model  calibration  is  to  obtain  a  mixed  path 
TD  calibration  model  which  is  consistent  with  available  model 
calibration  data  so  as  to  minimizes  the  difference  between  the 
measured  and  computed  TDs. 


3. 


CALIBRATION  DATA  ANALYSIS 


3 . 1  CALIBRATION  DATA 

The  U.S.  Coast  Guard-measured  data  provided  to  TASC 
for  model  calibration  include  TD  data  collected  at  land  (all¬ 
land  path)  sites  and  at  sea  (part  land  and  part  sea  water  path) 
sites.  The  land  data  set  includes  three  TDs/site:  TDw  (TDW) , 

TD  (TDX)  and  TD  (TDY) ,  collected  at  27  coastal  sites  distribu 
x  y 

ted  along  the  U.S.  West  Coast  and  shown  (by  triangles)  in  Fig. 
3.1-1.  The  land  site  locations  and  TD  measurements  are  as 
listed  in  Table  3.1-1.  The  sea  data  set  consists  of  two  TDs/ 
site  (TDX  and  TDY)  collected  at  23  sea  sites  located  in  the 
Southern  California  [between  Point  Arguello  (near  Santa 
Barbara)  and  San  Diego]  CCZ  as  shown  (by  circles)  in  Fig.  3.1-1 
The  sea  data  site  locations  and  TD  measurements  are  as  listed 
in  Table  3.1-2.  Based  on  information  from  discussions  with 
various  sources  and  engineering  judgement,  the  overall  quality 
of  the  TD  measurements  including  position  reference  errors  is 
assumed  to  be  between  0.1  -  0.2  psec. 


3.2  NONPARAMETRIC  DATA  ANLAYSIS 


The  purpose  of  nonparametric  analysis  of  the  calibra¬ 
tion  data  is  to  identify 


•  Potential  outliers  that  do  not  fit  the 
data  set 

•  Significant  functional  dependence(s)  in 
the  data  on  geophysical  propagation  pa¬ 
rameters 


R -48901 


Figure  3.1-1 


West  Coast  Loran-C  Chain  TD  Measurement 
Data  Site  Locations 
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NOTE:  Site  Coordinates  In  WGS-72. 

*No  measurement  data  available. 


TABLE  3.1-2 

SEA  DATA  SITE  LOCATIONS  AND  TIME  DIFFERENCE  MEASUREMENTS 

_ T-3550 


SITE 

COORDINATES 

TIME  DIFFERENCE  MEASUREMENTS  | 

IDENTIFICATION 

LATITUDE 

(NORTH) 

LONGITUDE 

(WEST) 

TDX 

TOY 

TASC 

USCG 

deg 

min 

sec 

deg 

min 

sec 

usee 

usee 

— 

2 

34 

43 

55*2 

■g 

44 

mrm 

27608.800 

41916.000 

1 

34 

42 

57*5 

47 

D 

27601.100 

41920.000 

7 

34 

26 

36*2 

120 

9 

27913.400 

41698.300 

3 

34 

26 

16*9 

120 

35 

27635.700 

41796.300 

9 

34 

23 

3.0 

1  19 

47 

32.6 

27976.400 

41592.400 

29 

6 

34 

16 

1.2 

120 

19 

32.6 

27680.300 

41699.000 

31 

5 

34 

11 

39ml 

119 

49 

14.1 

27965.400 

41546.000 

32 

4 

34 

5 

57.0 

119 

23 

1.4 

29036.200 

41410.500 

33 

8 

34 

5 

9.2 

120 

12 

48.8 

27696.000 

41615.700 

34 

13 

34 

2 

52*0 

119 

■1 

14.9 

28094.900 

41294.500 

35 

12 

34 

0 

48*3 

1  19 

o 

41.6 

26062.800 

41339.100 

36 

11 

33 

57 

ID 

D 

9.7 

26002.900 

41419.400 

37 

14 

33 

ID 

30 

15.0 

28165.400 

41122.500 

36 

10 

33 

1  19 

42 

23.6 

27976.600 

41450.900 

40 

IS 

33 

44*9 

1  18 

21 

29.9 

28174.400 

4! 036.000 

41 

16 

33 

44.5 

118 

10 

46.0 

28199.200 

40981.300 

42 

17 

33 

■ 

117 

47 

31.3 

28244.100 

40646. 100 

44 

19 

33 

117 

21 

46.  1 

28273.500 

406  75.300 

45 

16 

33 

I  1 6 

31 

15.0 

28120.400 

40975.600 

46 

20 

32 

55 

117 

17 

4.5 

28271.900 

40638.200 

47 

21 

32 

51 

^^9 

1  18 

37 

45.6 

28097.400 

40974.300 

48 

23 

32 

37 

46*6 

mm 

13 

1 6 1 0 

26262.6 00 

40601.000 

49 

22 

32 

37 

36.8 

m 

27.4 

28227.400 

40674.000 

NOTE:  No  TDW  measurement  data  are  available;  site  coordinates  in  WGS-72. 


•  Correlated  trends  between  land  and  sea 
data 

•  Likely  cause/effect  relationship  between 
data  and  geophysical  characteristics. 


For  each  TD  component  (i.e.,  TDW,  TDX  and  TDY),  the  land  and 
sea  subsets  were  analyzed  as  a  function  of 


•  Range  to  secondary  station 

•  Range  to  master  station 

•  Differential  range  between  secondary  and 

master  station 

•  Path  bearing  angle  at  master  station 

•  Path  bearing  angle  at  secondary  station. 

Figures  3.2-1  through  3.2-3  present  functional  dependence  plots 
of  TDW,  TDX  and  TDY  components,  respectively.  Each  figure 
shows  adjusted  TD  measurements  (Eq.  2.3-7)  as  a  function  of 
(a)  site  path  bearing  angle  (from  north)  measured  at  the  master 
station,  (b)  site  path  bearing  angle  (from  north)  measured  at 
the  secondary  station,  and  (c)  site  differential  range  (dif¬ 
ferential  primary  phase  delay)  between  secondary  and  master 
station,  as  defined  in  Fig.  A. 1-1  of  Appendix  A.  (Note,  the 
true  value  of  an  adjusted  TD  is  simply  the  difference  in  SFs 
for  the  secondary  and  master  station.  Since  sea  data  are  not 
available  for  TDW,  sea  data  plots  are  given  for  only  TDX  and 
TDY  components . ) 

In  addition  to  functional  dependence  trends,  examina¬ 
tions  of  these  figures  reveals  that  there  is  a  rather  strong 
correlation  between  land  and  sea  data  subsets  of  each  TD  com¬ 
ponent.  Figure  3.2-4  presents  a  composite  plot  of  all  three 
TD  components  for  the  ensemble  of  land  and  sea  data  as  a  func¬ 
tion  of  differential  range.  This  figure  is  presented  to  iden¬ 
tify  any  common  range-dependent  trend  (or  trends)  embodied  in 
all  of  the  three  TD  components. 

Analysis  of  measurement  data  based  on  plots  shown  in 
Figs.  3.2-1  through  3.2-4,  reveals  the  following: 
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ADJUSTED  TDW  (*i**c) 


3000  2000  1000  0  1000  2000  3000 

DIFFERENTIAL  RANGE  BETWEEN  SECONDARY 
AND  MASTER  STATION  (p**c) 


(c) 


Figure  3.2-1  Adjusted  TDW  (Land  Data)  as  a  Function  of 

(a)  Path-Bearing  Angle  at  Master  Station 

(b)  Path-Bearing  Angle  at  W  Secondary 
Station  and  (c)  Site  Differential  Range 
Between  W  Secondary  and  Master  Station 
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ADJUSTED  TDX  U**e, 


(a)  (b) 


DIFFERENTIAL  RANGE  BETWEEN  SECONDARY 
AND  MASTER  STATION  fci»e| 


LO 


Figure  3.2-2  Adjusted  TDX  (Land  and  Sea  Data)  as  a  Function 

of  (a)  Path  Bearing  Angle  at  Master  Station 
(b)  Path  Bearing  Angle  at  X  Secondary  Station 
and  (c)  Site  Differential  Range  Between  X 
Secondary  and  Master  Station 
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ADJUSTED  TOY  (mmc) 


BEARING  ANGLE  AT  SECONOARY  STATION  (iiaacl 


BEARING 


Figure  3. 


AT  MASTER  STATION  Ida?) 

(a) 


(b) 


ANO  MASTER  STATION  (vmcI 

(C) 


2-3  Adjusted  TDY  (Land  and  Sea  Data)  as  a  Function 
of  (a)  Path  Bearing  Angle  at  Master  Station 
(b)  Path  Bearing  Angle  at  Y  Secondary  Station 
and  (c)  Site  Differential  Range  Between  Y 
Secondary  and  Master  Station 


DIFFERENTIAL  RANGE  BETWEEN  SECONDARY  AND  MASTER  STATION  (pwc) 


Figure  3.2-4  Adjusted  TDW,  TDX  and  TDY  (both  Land 

and  Sea  Data)  Components  As  a  Function 
of  Site  Differential  Range  Between 
Secondary  and  Master  Station 


•  Dominant  and  very  similar  linear  range- 
dependent  trends  in  TDW  and  TDY  but  no 
identifiable  trend  in  TDX 

•  Highly  correlated  trend  between  land  and 
sea  data  as  a  function  of  both  range  and 
bearing  angle 

•  TDX  data  behavior  is  significantly  dif¬ 
ferent  than  that  of  TDW  or  TDY;  however, 
no  obvious  outliers  in  either  TDW,  TDX 
or  TDY  data 

•  TDX  data  identified  as  "A"  are  signifi¬ 
cantly  different  than  the  remainder  of 
TDX  data. 


3-9 


Further  examination  of  TDX  data  (although  not  essential  for 
the  semi-empirical  model  development)  indicated  Xhat  all  X 
station  radial  paths  in  the  TDX  data  between  bearing  angles  3 
and  150  deg  from  north  (identified  as  "A"  in  Figs.  3.2-2  (c)  and 
3.2-4)  exhibited  behavior  as  a  function  of  differential  range 
grossly  different  than  the  remainder  of  TDX  data.  This  sug¬ 
gested  the  possibility  that  a  terrain  with  propagation  proper¬ 
ties  different  from  those  of  the  remaining  chain  coverage  area 
may  exist  in  the  region  within  these  bearing  angles.  Indeed, 
the  location  of  the  San  Joaquin  Valley  (see  Fig.  3.1-1)  whose 
conductivity  is  higher  than  that  of  the  surrounding  area  by  an 
order  of  magnitude  (Ref.  10),  is  roughly  defined  by  this  region. 
Because  of  the  San  Joaquin  Valley's  orientation  relative  to 
the  X  station  and  the  shoreline,  all  X  station  radial  paths 
leading  to  the  CCZ  between  Los  Angeles  and  San  Diego  will  be 
significantly  impacted  by  the  presence  of  the  valley  and  will 
exhibit  signal  propagation  behavior  drastically  different  than 
the  rest  of  the  X  station  signal  coverage  area.  Note  that  the 
TDW  and  TDY  data  recorded  along  the  coastline  are  not  expected 
to  be  significantly  affected  by  the  high  conductivity  of  the 
valley.  This  is  because  the  propagation  path  segment  through 
the  valley  is  a  small  percentage  of  the  total  propagation  path 
for  the  M,  W,  and  Y  stations. 

3 . 3  SUMMARY 

Table  3.3-1  presents  a  summary  of  the  data  available 
for  calibrating  the  West  Coast  TD  grid  model.  Land  data  are 
distributed  along  the  West  Coast  from  Canada  to  San  Diego  while 
the  sea  data  are  concentrated  between  Point  Arguello  and  San 
Diego.  The  overall  quality  of  the  calibration  data  is  not 
known,  but  assumed  to  have  an  rms  accuracy  of  0.1  -0.2  psec. 
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TABLE  3.3-1 

WEST  COAST  MODEL  CALIBRATION  DATA  BASE  SUMMARY 


TYPE  OF 

NUMBER 

OF  DATA 
SITES 

NUMBER  OF 
MEASUREMENTS 

TOTAL  NUMBER 

OF  TD 

DATA  POINTS 

DATA 

TDW 

TDX 

TDY 

Land 

27 

25 

27 

24 

76 

Sea 

23 

* 

23 

23 

46 

Combined 
Land  and 
Sea 

50 

25 

50 

47 

122 

*No  data  available. 


The  data  are  considered  to  be  consistent  with  the 
expected  theoretical  behavior.  The  "apparently  anomalous" 
behavior  in  TDX  data  seem  to  have  been  caused  by  the  San 
Joaquin  valley  region  whose  conductivity  is  an  order  of  magni¬ 
tude  higher  than  the  surrounding  region. 

The  land  and  sea  TD  data  exhibit  similar  character¬ 
istics  as  a  function  of  both  station  range  and  bearing  angle 
at  the  station.  These  similar  characteristics  are  caused  by 
common  overland  path  segments  (i.e.,  conductivity )  in  land  and 
sea  data,  and  suggest  that  there  is  a  uniform  land/sea  water 
interface  vs  bearing  angle  effect  in  the  sea  data,  which  covers 
a  narrow  range  of  bearing  angles.  Additional  data  covering  a 
wider  range  of  bearing  angles  are  required  to  validate  the 
observed  uniformity  in  the  interface  effect.  The  observed 
data  characteristics  suggest  the  use  of  both  range  and  bearing 
angle  dependences  in  the  TD  model  structure.  Furthermore, 
especially  over  the  West  Coast  CCZ  sea  data  collection  region, 
these  characteristics  indicate  that  land  data  alone  may  be 
sufficient  to  calibrate  a  mixed  path  model. 


4. 


CALIBRATED  TD  MODEL 


4 . 1  INTRODUCTION 

The  purpose  of  this  chapter  is  to  present  an  accurate 
TD  grid  calibration  algorithm  for  the  West  Coast  Loran-C  chain 
CCZ.  Two  alternative  model  calibration  approaches,  shown  as  A 
and  B  in  Fig.  4.1-1,  are  considered.  Approach  A  is  designed 
to  assess  the  utility  of  using  only  land  based  data  for  CCZ 
model  calibration  as  compared  to  using  both  land  and  sea  based 
data  in  approach  B.  In  approach  A,  the  land  model  is  calibra¬ 
ted  while  the  sea  model  is  based  on  theory  (i.e.,  with  a  priori 
known  coefficients);  in  approach  B,  the  composite  land  and  sea 
model  is  calibrated  with  the  combined  land  and  sea  data.  Both 
approaches  apply  Millington's  empirical  method  to  combine  land 
and  sea  SFs  to  obtain  a  mixed  path  SF.  The  model  calibration 
procedure  for  both  approaches  is  detailed  in  the  next  section. 


4.2  TD  MODEL  CALIBRATION  EQUATIONS 

4.2.1  Adjusted  TD  Measurement  Equation 

It  is  convenient  to  model  the  transformed  form  of  the 
TD  measurements,  i.e.,  adjusted  TD  measurements,  given  by  Eq. 
2.3-8  and  repeated  below: 

ATD.  =  z.  =  (SF.  -  SFm)  +  vt  (4.2-1) 
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A-4MJI 


Figure  4.1-1  Alternative  Model  Calibration  Approaches 

where  SF .  (and  SF  )  are  the  SFs  associated  with  the  signal 
paths  from  i  (=w,  xory)  secondary  (and  master)  stations 
and  v^  is  the  total  measurement  error  associated  with  TD^. 

4.2.2  Sea  SF  Model 

The  structure  considered  for  the  sea  SF  model  is 

!®,i 

-tjF-  +  aQ  +  a^  T  (jsec,  if  10  <_  T  <_  540  (jsec 

(4.2-2) 

a  ,  «  . 

-~=—  +  a^  +  a1  T  Msec,  if  T  >  540  psec 

where  T  is  the  primary  phase  delay  (or  range);  a^  and  a^  (k  = 

-1,  0  and  1)  are  sea  coefficients  considered  as  known  (no  un¬ 
certainty)  for  approach  A  and  unknown  (uncertain)  for  approach 
B.  The  sea  model  coefficient  values  for  approach  A  are  (Ref.  9) 
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a_2  =  2.741 

a  =  -0.0114 
o 

ax  =  0.0003277  (4.2-3) 

al±  =  129.043 

a'  =  -0.408 
o 

aj  =  0.0006458 
4.2.3  Land  SF  Model 

A  general  polynomial  model  structure  for  the  land  SF 
associated  with  a  station  is  given  by  Eq.  2.2-4.  By  special¬ 
izing  this  model  structure  to  the  West  Coast  chain  CCZ  service 
area  and  chain  topography,  and  incorporating  the  results  of 
calibration  data  analysis,  the  following  two  candidate  forms 
for  the  land  SF  model  were  considered: 


“Localized”  Range/Bearing  (LRB)  Model 


SFj  =  Aq  +  [Ax  +  E1  fx(Px)l  T.  Msec  (4.2-4) 


"Generalized”  Range/Bearing  (GRB)  Model 


SFj  =  Ao  +  A1  Tj  +  [Cj£  sin  £pj 


+  cos  £0 . J  Msec 


(4.2-5) 


where  A  ,  A,,  B, ,  C.0  and  D.«  are  uncertain  model  coefficients. 

Oil  JX  JX 

p.  is  the  path  bearing  angle  at  the  j  (=  w,  x,  y  or  m)  sta- 
J  «-k 


tion  and  Tj  is  the  path  range  to  the  j 


station;  function 


fx(Px)  is  zero  for  all  chain  stations  except  for  the  X  station. 
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For  the  X  station,  f  is  zero  unless  the  X  station  radial  (sig 
nal  path)  passes  through  the  San  Joaquin  Valley,  then  it  is 
unity . 


The  LRB  model  is  purposely  kept  as  simple  (fewest  un¬ 
certain  model  coefficients)  as  possible  yet  designed  to  embody 
distortions  (warpages)  to  the  X  station  SF  caused  by  the  San 
Joaquin  valley.  The  GRB  model,  on  the  other  hand,  includes 
bearing  angle  dependences  for  all  four  chain  stations  instead 
of  the  X  station  alone,  as  is  the  case  in  the  LRB  model.  Con¬ 
sequently,  the  GRB  model  is  relatively  more  complex  and  is 
expected  to  exhibit  performance  superior  to  the  LRB  model. 

Note  that  a  calibrated  model  is  expected  to  be  accurate  and 
applicable  only  over  the  extent  of  ranges  and  bearing  angles 
embodied  in  the  calibration  data.  Hence,  outside  the  region 
covered  by  the  calibration  data,  the  model  may  not  be  as  accu¬ 
rate  as  within  the  data  coverage  region. 

4. 2. 4  Mixed  Path  SF  Computations 

The  SF  over  a  mixed  path  is  computed  using  Millington's 
empirical  equations,  Eqs.  2.2-1  through  2.2-3.  In  these  equa¬ 
tions,  all  terms  except  SF(o^,  Tg)  in  Eq.  2.2-3  can  be  computed 
with  land  and  sea  SF  models  described  in  Sections  4.2.2  and 
4.2.3.  The  SF(o^»  Tg)  is  the  SF  of  a  fictitious  land  path  of 
length  Tg  (the  sea  segment  of  the  actual  mixed  path)  which  is 
usually  much  smaller  than  any  of  the  actual  land  segment  path 
lengths  embodied  in  the  West  Coast  land  calibration  data  base. 
Consequently,  the  land  SF  model  to  be  developed  with  the  cali¬ 
bration  data  base  cannot  be  used  to  compute  the  term  SF(a^»  Tg) 
in  Eq.  2.2-3.  Therefore,  the  following  theoretical  polynomial 
land  SF  model  is  used: 

SF(ol,  Tg)  =  +  0.439  +  0.00245  Tg  psec  (4.2-6) 
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where  Tg  is  in  psec.  This  polynomial  model  was  derived  by 
fitting  it  to  the  homogeneous/smooth  earth  theoretical  predic¬ 
tions  (Ref.  4).  The  coefficients  of  this  model  correspond  to 
an  average  ground  conductivity  of  0.003  mhos/m,  which  is  the 
estimated  average  conductivity  of  the  West  Coast  chain  coverage 
area  based  on  the  nonparametric  data  analysis  results  (Chapter 
3)  and  homogeneous/smooth  earth  theory  (Ref.  4).  Mixed  path 
SFs  derived  from  the  combination  of  land  and  sea  SF  models,  as 
per  Eqs.  2.2-1  through  2.2-3,  are  then  incorporated  in  the 
adjusted  TD  measurement  equation,  Eq.  4.2-1,  which  is  cali¬ 
brated  with  the  measurement  data  as  described  in  the  next  sec¬ 
tion. 


4.3  MODEL  CALIBRATION 

4.3.1  Calibration  Procedure 

An  overview  of  the  West  Coast  TD  grid  calibration 
procedure  is  illustrated  in  Fig.  4.3-1.  The  first  step  in  the 
model  calibration  procedure  is  to  hypothesize  candidate  model 
structures  for  the  land  and  sea  SF  models.  In  approach  A,  the 
sea  model  is  known  (theoretical)  and  therefore  only  the  land 
model  is  hypothesized  as  opposed  to  approach  B  where  both  land 
and  sea  models  are  postulated. 

Next,  mixed  path  SFs  are  substituted  to  form  adjusted 
TD  measurements  (Eq.  4.2-1)  which  can  be  symbolically  written 
in  matrix  form  (Ref.  3)  as 

z  =  H  x  +  v  (4.3-1) 
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W  —0*0 


Figure  4.3-1  TD  Grid  Calibration  Procedure 


where  the  measurement  vector  is 


z 


ATD 

w 

ATD 

x 

ATD 

y 


(4.3-2) 


The  observation  matrix,  H,  is  a  function  of  range  and  bearing 
angle  of  station  paths  associated  with  the  TD.  The  state  vec¬ 
tor,  x,  is  a  vector  of\uncertain  model  coefficients  which  are 
to  be  estimated  from  daia.  The  TD  measurement  error  vector  is 
v  given  by 


w 


v 


V 


y 


(4.3-3) 
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where  TD  component  (i.e.,  TDW,  TDX  or  TDY),  and  hence  ATD 
component,  measurement  errors  are  assumed  to  be  random  with 
zero  mean. 

The  next  step  in  the  calibration  procedure  is  to  pro¬ 
vide  a  priori  information  on  the  initial  estimates  and  uncer¬ 
tainties  of  model  coefficients,  and  TD  measurement  error  sta¬ 
tistics  to  the  coefficient  estimation  algorithm.  The  Kalman 
filter  (Ref.  3)  provides  a  convenient  method  to  estimate  the 
state  vector  (coefficients)  with  a  zero  mean  using  the  adjusted 
TD  measurements.  The  a  priori  information  was  developed  by  a 
combination  of  data  anlaysis  results  and  expected  theoretical 
behavior  of  Loran-C  signals  over  land  and  sea  water  paths.  In 
particular,  the  sea  water  SF  model  coefficients  were  constrained 
to  reasonable  theoretical  limits  to  warrant  the  calibrated 
model  useful  beyond  the  CCZ  region  (where  no  data  were  avail¬ 
able  for  model  calibration). 

Initial  attempts  to  calibrate  candidate  West  Coast  TD 
models  consistently  yielded  TD  residuals  with  an  rms  level  of 
about  0.4  psec.  Therefore,  TD  component  measurement  error  in 
the  calibration  procedure  was  assumed  to  be  0.4  psec.  Note, 

0.4  psec  error  includes  receiver  measurement  error  of  0.1  - 
0.2  psec  (Chapter  3),  data  site  position  location  reference 
errors,  unmodeled  TD  warpage  conditions  and  emission  delay 
variations  from  site-to-site . 

More  than  20  candidate  GRB  model  structures  (Eq.  4.2-5) 
with  varying  numbers  of  harmonic  terms  as  well  as  the  single 
LRB  model  structure  (Eq.  4.2-4)  were  considered  for  the  land 
SF  model.  The  sea  SF  model  was  always  chosen  to  be  3-term, 
range -dependent  model  given  by  Eq.  4.2-2. 

Each  candidate  land  SF  model  in  combination  with  the 
theoretical  sea  SF  model  (approach  A),  or  semi-empirical  sea  SF 
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model  (approach  B),  was  calibrated  with  land  data  (approach  A) 
or  combined  land  and  sea  data  (approach  B).  For  each  candidate 
combination,  an  associated  candidate  observation  matrix  and  a 
state  vector  were  developed. 

Performance  of  each  calibrated  candidate  combined 
land  and  sea  model  was  evaluated  in  terms  of  the  statistical 
reasonableness  of  the  calibrated  model  fit  to  the  data.  Sta- 
tisitacal  reasonableness  was  quantified  as  indicated  in  Fig. 
4.3-2  in  terms  of  the  standard  deviation  of  individual  TD  com¬ 
ponent  site  residuals  and  the  chi-square  test  which  evaluates 
the  fit  of  the  three-  (or  two-)  dimensional  TD  site  residual 
(i.e.,  all  TD  components  together  at  a  site)  to  the  data. 

4.3.2  Calibrated  Models 

Of  the  candidate  GRB  land  SF  models  considered,  the 
TD  model  with  the  following  GRB  functional  form  yielded  the 
best  performance: 


sfj 


A.  +  A,  T. 
O  1  J 


+  !?1  [CJ*  Si"  +  V  C°S 


(4.3-4) 


where 


C 

D 

C 

D 


j  =  w,  x, 

wl  =  Cxl  = 
wl  “  Dxl 


y 

c 

D 


or  m  station 


yi  ' 
yi  = 


o 

o 


Table  4.3-1  describes  the  number  of  uncertain  coefficients  in¬ 
cluded  in  the  "best"  GRB  model  and  the  LRB  model  under  each  of 
the  two  calibration  approaches  considered. 
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TABLE  4.3-1 

NUMBER  OF  COEFFICIENTS  IN  LRB  AND  GRB  TD  MODELS 


Note,  in  approach  A  where  only  land  data  are  used  to 
calibrate  the  land  model,  one  bias  state  (coefficient)  per 
land  TD  component  is  required  in  the  TD  model  (see  Table  4.3-1) 
to  account  for  a  possible  constant  bias  (shift)  in  the  secondary 
station  emission  delay  and  the  unobservable  biases  in  the  land 
SF  models.  However,  in  approach  B  where  both  land  and  sea 
data  are  used  to  calibrate  the  composite  land  and  sea  model, 
an  additional  bias  state  per  sea  TD  component  is  included  as 
shown  in  Table  4.3-1)  to  result  in  zero-mean  residuals.  If 
sea  TD  bias  states  are  not  included  in  the  model,  the  mean  TD 
residual  is  non-zero,  and  furthermore  the  rms  TD  residual  is 
significantly  larger  than  that  obtained  with  sea  TD  bias  states. 
Thus,  a  total  of  three  land  TD  bias  states  (for  TDW,  TDX  and 
TDY)  are  included  in  the  TD  model  calibrated  in  approach  A, 
while  three  land  and  two  sea  TD  bias  states  (for  TDX  and  TDY) 
are  used  in  the  TD  model  calibrated  in  approach  B.  As  expected, 
the  magnitude  of  the  sea  bias  states  in  the  model  calibrated 
in  approach  B  were  roughly  the  same  as  the  corresponding  means 
in  the  sea  TD  residuals  obtained  in  approach  A  where  sea  TDs 
are  computed  using  a  calibrated  land  model  and  the  theoretical 
sea  model,  Eqs.  4.2-2  and  4.2-3.  There  is  not  sufficient  data 
to  identify  the  likely  sources  of  the  observed  sea  TD  biases. 

The  land/sea  water  interface  "phase  recovery"  effect  (Ref.  7) 
may  be  responsible  for  part  of  the  observed  sea  TD  biases. 


4.4  CALIBRATED  MODEL  PERFORMANCE 

4.4.1  Calibration  Data  Base 

Table  4.4-1  summarizes  the  rms  TD  residuals  for  the 
LRB  and  GRB  TD  models  calibrated  with  (1)  land  data  alone  (ap¬ 
proach  A)  and  (2)  combined  land  and  sea  data  (approach  B).  In 
this  table,  the  rms  value  of  the  residuals  at  the  46  sea 
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TABLE  4.4-1 

PERFORMANCE  COMPARISON  OF  LRB  AND  GRB  MODELS 


T-3566 


CALIBRATION 

APPROACH 

MODEL 

NUMBER  OF 
CALIBRATION 

DATA  POINTS 

RMS  TD  RESIDUAL  -  pSEC 

OVEK  SEA 
DATA  POINTS 
(46) 

OVER  COMBINED 
LAND  AND  SEA 
DATA  POINTS 
(122) 

LAND 

SEA 

TOTAL 

A 

ibrated 
d  Hodel 
eoretical 
Model) 

LRB 

GRB 

76 

76 

- 

76 

76 

0.703 

0.769 

0.636 

0.570 

Cal 

Lan 

(Th 

Sea 

B 

LRB 

GRB 

76 

76 

46 

46 

122 

122 

■ 

0.521 

0.390 

Calibrated 

Composite 

Land  and 

Sea  Hodel 

calibration  TD  data  points  (and  also  for  over  the  entire  set 
of  122  land  and  sea  calibration  TD  data  points)  are  presented. 

Comparison  of  sea  residual  statistics  (Table  4.4-1) 
obtained  with  the  LRB  and  GRB  models  shows  similar  performance 
for  both  models  in  either  approach.  However,  calibration  ap¬ 
proach  B  yields  a  factor  of  two  improvement  in  the  sea  rms 
residual  over  those  obtained  in  approach  A.  Thus,  from  con¬ 
siderations  of  sea  residuals  alone,  approach  B  is  preferred 
over  approach  A.  Further  comparison  of  rms  residuals  obtained 
with  the  LRB  and  GRB  models  (Table  4.4-1)  over  the  entire  set 
of  land  and  sea  calibration  data  points  indicates  that  the  GRB 
model  yields  superior  performance.  Therefore  the  GRB  model 
calibrated  with  land  and  sea  data  was  selected  as  the  ''best11 
performance  TD  model  for  calibrating  the  West  Coast  Loran-C 
chain.  The  residual  statistics  obtained  with  the  GRB  model 
for  the  indicated  TD  components  and  data  sets  are  summarized 
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in  Table  4.4-2.  The  details  of  the  calibrated  GRB  model  and 
LRB  model  algorithms  are  given  in  Appendix  A. 


TABLE  4.4-2 

TD  RESIDUAL  STATISTICS  OF  GRB  MODEL 
OVER  CALIBRATION  DATA  BASE  (APPROACH  B) 


CALIBRATION 
DATA  BASE 

RMS  TD  RESIDUAL  -  psec 

TDW 

TDX 

TDY 

COMBINED  TD 
COMPONENTS 

Land 

0.308 

0.501 

0.382 

0.408 

Sea 

* 

0.388 

0.306 

0.350 

Combined 

Land  and  Sea 

0.308 

0.457 

0.347 

0.390 

*No  calibration  sea  data  available. 


Figure  4.4-1  presents  the  TD  residuals  (solid  curve) 
at  each  land  and  sea  data  site  obtained  with  the  calibrated 
West  Coast  GRB  model.  The  data  collection  sites  are  arranged 
in  order  from  north  to  south.  For  comparison,  Fig.  4.4-1  also 
shows  the  calibration  data  (dotted  line)  and  the  adjusted  TD 
measurements.  Both  data  and  residuals  have  breaks  (or  gaps) 
at  sites  where  no  measurement  data  are  available.  (Note,  there 
are  no  sea  data  for  the  TDW  component.) 

4.4.2  Validation  Data  Base 

The  calibrated  model  was  also  evaluated  at  the  25  sea 
sites  listed  in  Table  4.4-3,  the  coordinates  of  which  were 
specified  by  the  U.S.  Coast  Guard.  The  computed  TDs  at  the 
validation  sites,  as  listed  in  Table  4.4-3,  were  forwarded  to 
the  U.S.  Coast  Guard  for  comparison  with  TD  measurements.  The 
rms  TD  residual  over  the  ensemble  of  all  TDs  (i.e.,  25-TDX 
and  25-TDY)  computed  by  the  U.S.  Coast  Guard  is  0.42  psec,  a 
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(a)  TPW  Residual 


(b)  TDX  Residual 


0  1*  »  M  «0  M 

SITf  NUMtfft 

Non:  Ink  TO  Mownom  m  *  on  •  *"#*e*wd  by  •  pp  m  t*»  otot 


(c)  TDY  Residual 


Figure  4.4-1 


TD  Residuals  for  GRB  Model  Calibrated  with 
Combined  Land  and  Sea  Data  (Approach  B) 
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factor  of  four  improvement  over  the  original  U.S.  Coast  Guard 
TD  grid  charting  procedures  (Ref.  10).  Note,  the  rms  residual 
over  the  validation  data  points  (0.42  psec)  is  approximately 
the  same  value  as  that  obtained  over  the  calibration  data 
points  (0.35  psec). 

4 . 5  SUMMARY 

Land  data  alone  (approach  A)  can  be  used  to  calibrate 
the  West  Coast  CCZ  grid  to  an  rms  accuracy  of  approximately 
0.8  psec.  Use  of  sea  data  in  addition  to  land  data  (approach 
B)  provides  a  factor  of  two  improvement  in  the  CCZ  grid  accuracy. 
The  GRB  model  calibrated  with  combined  land  and  sea  data  yields 
the  best  overall  rms  TD  residual  performance  and  was  therefore 
selected  as  the  West  Coast  calibrated  TD  grid  model.  The  cali¬ 
brated  GRB  model  has  an  rms  TD  error  of  0.35  psec  over  sea 
calibration  data  points  which  would  result  in  an  rms  position 
error  of  340  m  in  the  West  Coast  CCZ  if  TDX  and  TDY  LOPs  are 
utilized  for  the  position  fix  by  a  user  with  a  "perfect"  re¬ 
ceiver.  The  rms  TD  error  of  the  calibrated  model  over  the 
U.S.  Coast  Guard  sea  validation  data  points  (not  used  in  model 
calibration)  is  0.42  psec  resulting  in  a  factor  of  four  im¬ 
provement  over  the  original  charting  procedures.  This  demon¬ 
strates  that  semi-empirical  grid  calibration  techniques  are 
effective  for  calibrating  an  accurate  Loran-C  grid  for  the 


Coastal  Confluence  Zone. 


5. 


MODEL  ACCURACY  SENSITIVITY  ANALYSIS 


5 . 1  INTRODUCTION 

The  sensitivity  analysis  presented  herein  assesses 
the  operational  practicality  of  adopting  semi-empirical  tech¬ 
niques  (or  models)  as  a  TD  grid  calibration  tool  for  Loran-C 
chain  CCZ  regions.  Key  operational  issues  are  the  quantity 
and  distribution  of  grid  calibration  data  required  to  achieve 
a  desired  grid  accuracy.  This  section  examines  the  accuracy 
of  the  West  Coast  semi-empirical  GRB  TD  model  grid  (developed 
in  Chapter  4)  in  terms  of  the  quantity  and  distribution  of 
data  used  for  calibrating  the  model.  Based  on  West  Coast  sen¬ 
sitivity  analysis  results,  calibration  data  collection  guide¬ 
lines  are  formulated  to  aid  in  the  design  of  future  semi- 
empirical  grid  calibration  efforts. 


5.2  ANALYSIS 

5.2.1  Sensitivity  and  Evaluation  Data  Bases 

*  Two  mutually  exclusive  data  bases  (Table  5.2-1)  are 
formed,  for  the  sensitivity  analysis  studies,  out  of  the  avail¬ 
able  TD  measurement  data  used  in  Chapter  4  to  calibrate  the 
West  Coast  TD  grid  model.  The  two  data  bases  are  refered  to 
as  the  sensitivity  data  base  and  the  evaluation  data  base. 

The  sensitivity  data  base  consists  of  both  land  and  sea  data 
sites  (distributed  from  Canada  to  San  Diego),  subsets  of 
which  are  used  to  calibrate  the  sensitivity  analysis  model. 

The  accuracy  of  each  calibrated  model  is  assessed  with  the 


TABLE  5.2-1 

SENSITIVITY  AND  EVALUATION  DATA  BASES  SUMMARY 


DATA  BASE 

NUMBER  OF 

DATA  SITES 

NUMBER  OF 

TD  DATA  POINTS 

LAND 

SEA 

TOTAL 

LAND 

SEA 

TOTAL 

Sensitivity 

27 

12 

39 

76 

24 

100 

Evaluation 

- 

11 

11 

- 

22 

22 

Combined 

27 

23 

50 

76 

46 

122 

evaluation  data  base  which  includes  only  sea  sites  distributed 
in  the  Southern  California  CCZ  (between  Point  Arguello  and 
San  Diego) . 


5.2.2  Approach 


A  number  of  subsets  of  the  sensitivity  data  base  are 
used  to  calibrate  the  sensitivity  model.  These  subsets  in¬ 
cluded 


•  Uniform  distributions  of  combined  land 
and  sea  calibration  data  sites  as  shown 
in  Fig.  5.2-1 

•  Clusters  of  land  calibration  sites  (lo¬ 
cated  either  north  or  south  of  San  Fran¬ 
cisco)  with  uniform  distribution  of  sea 
sites  located  in  Southern  California,  as 
shown  in  Fig.  5.2-2. 


Each  of  the  calibration  data  sets  (subsets  of  the 
sensitivity  data  base)  are  used  (one  at  a  time)  to  calibrate 
the  sensitivity  model  by  the  calibration  procedure  described 
in  Section  4.  The  calibrated  model  is  then  used  to  compute 
the  TDs  and  resulting  TD  residuals  at  all  sites  in  the  combined 
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.2-1  Illustration  of  Clustered  Land  Calibration 
Sites  and  Uniformly  Spaced  Sea  Calibration 
Sites 


Figure  5.2-2 


Illustration  of  Uniformly  Distributed 
Combined  Land  and  Sea  Calibration  Data  Sites 


data  base.  The  rms  of  the  TD  residuals  of  each  calibrated 
model  is  then  computed  over  (1)  the  "evaluation  sites"  (i.e., 
ensemble  of  all  TD  components  at  all  sites  in  the  evaluation 
data  base)  and  (2)  "all  sites"  (i.e.,  ensemble  of  all  TD  com¬ 
ponents  at  all  sites  in  the  combined  data  base). 

5.2.3  Sensitivity  Analysis  Results 

The  rms  TD  residuals  obtained  with  the  GRB  TD  model 
calibrated  with  subsets  of  the  sensitivity  data  base  are  shown 
in  Fig.  5.2-3  as  a  function  of  quantity  and  distribution  of 
the  calibration  data,  where 


•  The  quantity  of  calibration  data  is  ex¬ 
pressed  as  a  percentage  of  TD  data  points 
in  the  combined  data  base  (which  includes 
data  points  in  both  sensitivity  and  eval¬ 
uation  data  bases) 

•  The  data  distribution  is  keyed  as  a  bar 
(representing  uniform  data  site  distri¬ 
bution)  or  a  triangle  (denoting  clusters 
of  calibration  data  sites) 

•  The  length  of  a  bar  shows  the  spread  of 
the  computed  rms  TD  residuals  obtained 
for  models  calibrated  with  several  uni¬ 
formly  distributed  subsets  of  the  senfi- 
tivity  data  base,  each  containing  approxi¬ 
mately  the  same  number  of  data  points 

•  Adjacent  solid  and  open  areas,  bars  or 
triangles,  are  the  corresponding  rms  TD 
residuals  over  all  sites  and  evaluation 
sites,  respectively. 


The  rms  residuals  shown  for  100  percent  of  the  data  are  the 
residuals  of  the  model  calibrated  with  all  the  available  data 
in  Chapter  4  (i.e.,  -including  evaluation  data)  and  are  shown 
for  comparison  purposes.  Brief  explanations  of  the  observed 
calibrated  model  accuracy  behavior  for  both  clustered  and 
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Figure  5.2-3  Model  Accuracy  Sensitivity  to  Quantity 

and  Distribution  of  Calibration  Data 


uniformly  distributed  calibration  data  conditions  are  presented 
in  the  following  sections. 

5. 2. A  Clustered  Sets  of  Calibration  Data 


Two  specific  examples  of  clustered  data  sets  are  shown 
in  Fig.  5.2-3.  Calibration  data  set  A  (see  Fig.  5.2-1)  includes 
all  land  data  south  of  San  Francisco  and  all  sea  data  available 
in  the  sensitivity  data  base  (also  south  of  San  Francisco). 

No  calibration  data  north  of  San  Francisco  is  contained  in 
data  set  A.  Consequently,  the  model  calibrated  with  data  set 
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A  performs  poorly  over  data  sites  north  of  San  Francisco  as 
manifested  by  the  high  rms  residual  computed  over  all  sites. 

On  the  other  hand,  the  calibrated  model  is  expected  to  exhibit 
excellent  performance  in  the  region  south  of  San  Francisco 
where  all  of  the  evaluation  sites  are  located.  This  perform¬ 
ance  is  manifested  by  the  low  residual  value  computed  over  the 
evaluation  sites. 


Calibration  data  set  B  includes  land  data  north  of 
San  Francisco  and  all  sea  data  (south  of  San  Francisco)  -- 
thus,  the  sensitivity  data  base  includes  data  distributed  (al¬ 
though  not  uniformly)  over  the  entire  West  Coast.  Consequently, 
the  model  calibrated  with  data  set  B  is  expected  to  have  a 
relatively  better  residual  performance,  especially  over  the 
data  sites  located  north  of  San  Francisco,  than  the  model  cali¬ 
brated  with  A.  This  is  manifested  in  Fig.  5.2-3  by  a  smaller 
rms  residual  over  all  sites  for  case  B.  Note,  evaluation  site 


residuals  for  both  cases  are  comparable,  since  both  subsets  A 
and  B  include  calibration  data  over  the  region  covered  by  the 
evaluation  sites.  From  the  comparison  of  calibrated  model 
performance  for  the  two  sets,  it  is  concluded  that  the  calibra¬ 


tion  data  set  must  be  representative  of  the  region  to  be  cali¬ 


brated. 


5-2.5  Uniformly  Distributed  Calibration  Data  Sets 


A  number  of  uniformly  distributed  calibration  data 
sets  with  varying  data  density  were  analyzed.  As  an  example, 
consider  the  computed  spread  of  residuals,  labeled  as  C  in 
Fig.  5.2-3.  Adjoining  residual  bars  labeled  all  sites  and 
evaluation  sites  correspond  to  the  use  of  18  calibration  data 
points,  roughly  13  percent  of  the  combined  data.  These  bars 
("C")  depict  the  spread  of  computed  rms  residuals  for  four 
different  uniformly  distributed  calibration  data  subsets  formed 
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from  the  sensitivity  data  base  by  retaining  every  sixth  data 
site.  For  this  case,  two  of  the  four  subsets  of  the  sensitiv¬ 
ity  data  base  considered  are  labeled  1  and  2  in  Fig.  5.2-1. 
Subsets  1  and  2  are  similar  except  subset  2  does  not  span  the 
northern  tip  of  the  U.S.  West  Coast  area  (identified  as  "U"  in 
Fig.  5.2-2)  and  covered  by  sites  1  through  4.  Thus,  the  model 
calibrated  with  subset  2  extrapolates  over  "U"  while  the  model 
calibrated  with  subset  1  interpolates  over  "U" .  Therefore, 
the  all  site  residual  performance  of  the  model  calibrated  with 
subset  2  is  inferior  to  the  performance  using  subset  1 ,  as 
manifested  by  the  highest  all  site  residual  value.  There  is 
very  little  spread,  as  expected,  in  the  evaluation  site  resid¬ 
uals  because  all  four  calibration  data  subsets  cover  the  region 
of  evaluation  sites. 

As  expected,  the  spread  in  rms  residual  values  (Fig. 
5.2-3)  associated  with  both  all  sites  and  with  the  evaluation 
sites  decreases  with  increasing  density  of  the  data  in  a  uni¬ 
formly  distributed  calibration  data  set.  No  significant  im¬ 
provement  in  the  West  Coast  calibrated  grid  rms  residual  per¬ 
formance  is  observed  with  the  use  of  more  than  40  percent  of 
the  available  data  for  model  calibration  data. 


5 . 3  SUMMARY 


Based  on  the  West  Coast  model  accuracy  sensitivity 
analysis  results,  it  is  concluded  that 


•  A  uniform  distribution  of  50  percent  of 
the  available  measurement  data  provides 
a  CCZ  TD  grid  with  an  rms  accuracy  of 
0.4  Msec  as  compared  to  0.35  psec  with 
100  percent  of  data 


•  Data  collection  sites  should  be  selected 
to  provide  a  relatively  uniform  distribu¬ 
tion  along  the  coast  with  an  average 
spacing  of  100-200  km  over  the  region  of 
interest  (as  compared  to  50-100  km  for 
land  sites  and  ~20  km  for  sea  sites  in 
the  available  measurement  data) 

•  Additional  data  sites  should  be  concen¬ 
trated  in  regions  receiving  signal  paths 
through  known  or  suspected  anomalous 
propagation  region(s)  (e.g. ,  San  Joaquin 
Valley) . 

•  Combination  of  land  and  sea  data  yields 
a  higher  accuracy  calibrated  grid  than 
possible  with  land  data  alone. 


Although  the  issue  of  utilizing  land  vs  sea  data  has  not  been 
fully  investigated  due  to  limited  quantity  and  spatial  (coastal) 
coverage  provided  by  the  available  sea  calibration  data,  pre¬ 
liminary  results  indicate  the  following: 


•  Either  land  or  sea  data  may  be  used  to 
calibrate  a  CCZ  grid 

•  Some  sea  data  are  always  desired  to  iden¬ 
tify  land/sea  interface  effects 

•  Sea  data  collected  over  a  wider  coastal 
region  will  help  to  identify  potential 
source(s)  of  the  sea  bias  seen  in  the 
present  study 

•  Both  near  and  far  from  shore  sea  data 
will  provide  greater  observability  to 
sea  model  parameters  (coefficients) 

•  Inclusion  of  land  data  greatly  reduces 
the  required  density  of  sea  calibration 
data  as  the  dominant  bearing  dependence 
effects  are  easily  observable  in  the 
land  data. 


It  is  recommended  that  additional  sea  data  of  the  type  described 
above  be  used  to  verify  the  above  preliminary  findings. 


6. 


CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  CONCLUSIONS 


The  utility  of  semi-empirical  techniques  to  accurately 
calibrate  Loran-C  grids  in  the  CCZ  has  been  demonstrated  by 
applying  these  techniques  to  the  West  Coast  Loran-C  CCZ.  In 
this  region,  current  U.S.  Coast  Guard  prediction  procedures 
have  been  reported  to  result  in  significant  charting  errors 
(Ref.  10).  The  TASC-deve loped  algorithm  exhibits  the  follow¬ 
ing  characteristics: 


•  *  RMS  TD  error  of  0.42  psec  over  the  sea 

TD  measurement  data  points  not  used  in 
model  calibration  and  0.35  psec  over  the 
data  points  used  in  model  calibration 

•  RMS  position  error  of  less  than  400  m  in 
the  Southern  California  CCZ  --  a  factor 
of  four  improvement  over  the  original 
charting  procedures  (Ref.  10T 

•  Reasonably  accurate  TDs  beyond  the  CCZ 
(where  measurement  data  were  not  avail¬ 
able  for  model  calibration) 

•  Computationally  simple  (can  be  implemen¬ 
ted  on  a  hand-held  electronic  calculator 
similar  to  the  HP-67) 

•  Cost  effective  since  much  less  calibra¬ 
tion  data  are  required  than  for  other 
known  calibration  procedures 

•  Can  be  easily  extended  to  include  data 
from  other  coverage  regions  as  it  be¬ 
comes  available. 


In  addition  the  results  of  this  study  show: 


•  Land  data  alone  can  be  used  to  calibrate 
the  West  Coast  CCZ  grid  with  an  rms  TD 
error  of  approximately  0.8  psec 

•  Inclusion  of  sea  calibration  data  pro¬ 
duces  a  factor  of  two  improvement  in  the 
calibrated  grid  accuracy 

o  A  uniform  distribution  of  50  percent  of 
the  available  data  (average  spacing  of 
100-200  km)  provides  an  rms  CCZ  grid 
accuracy  of  0.40  psec  as  compared  to 
0.35  psec  acheived  with  100%  of  the 
available  data. 


In  summary,  the  semi-empirical  TD  grid  calibration  techniques 
have  been  shown  to  be  both  effective  and  efficient  for  develop¬ 
ing  accurate  CCZ  grid. 


6 . 2  RECOMMENDATIONS 

It  is  recommended  that  semi-empirical  grid  calibration 
techniques  using  both  land  and  sea  calibration  data  be  applied 
to  other  Loran-C  chains  and  regions  of  the  Coastal  Confluence 
Zone  to  develop  accurate  TD  grids.  In  particular,  it  is  recom¬ 
mended  that  this  technique  be  applied  to  develop  accurate  TD 
grids  for: 

•  Great  Lakes 

•  East  Coast 

•  Gulf  of  Mexico. 

Also,  it  is  recommended  that  a  grid  calibration  data  collection 
and  management  program  plan  be  established  to 
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•  Design  data  collection  requirements  for 
future  semi-empirical  grid  calibration 
of  Loran-C  chains 

•  Develop  procedures  and  methods  for  col¬ 
lecting  the  most  useful  data  and  only 
necessary  data 

•  Provide  a  capability  for  on-line  inter¬ 
action  with  the  data  collection  team  to 
identify  and  verify  "abnormal"  data  be¬ 
havior 

•  Manage  the  collected  data  so  as  to  pro¬ 
vide  reliable  and  efficient  computer 
access  to  any  set  or  subset  of  raw  or 
processed  data 

•  Develop  data  handling  and  analysis  soft¬ 
ware. 


The  data  collection  and  management  program  plan  outlined  above 
wil  provide  the  U.S.  Coast  Guard  with  a  cost-effective  technique 
for  semi-empirical  Loran-C  grid  calibration  and  chart  validation. 


APPENDIX  A 


CALIBRATED  TD  GRID  ALGORITHMS 

A.  1  INTRODUCTION 

This  appendix  presents  the  TD  grid  algorithms  for  the 
two  "finalist"  calibration  models,  identified  as  the  GRB  model 
and  LRB  model.  The  GRB  model  has  been  selected  as  the  West 
Coast  TD  grid  calibration  model.  Table  A. 1-1  gives  a  computa¬ 
tional  guide  to  equations  and  tables  that  are  required  to  com¬ 
pute  TDs,  which  are  presented  in  this  appendix.  The  TD  is  the 
difference  between  the  times-of-arrival  of  signals  from  the  i*"*1 
secondary  (w,  x  or  y)  and  master  (m)  stations  at  a  user  as  il¬ 
lustrated  in  Fig.  A. 1-1  and  expressed  by  the  following  equation 

TD.  =  (T.-T)  +  (SF.-SF)  +  ED.  +  b,  (A. 1-1) 

x  i  m  i  m  lx 

where 

n  R. 

T.  =  Msec  (A. 1-2) 

n  R 

Tm  =  MSec  (A. 1-3) 

R.  =  i£k  secondary  station- to-user  great- 
circle  path  length 

R  =  master  station-to-user  great-circle  path 
m  length 
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TABLE  A. 1-1 
TD  COMPUTATION  GUIDE 


TD  Equation:  TD.  =  T.  -  T  +  SF .  -  SF  +  ED.  +  b. 

limimii 

TD  EQUATION  TERMS 

EQUATION 

NUMBER 

TABLE 

NUMBER 

Ti 

A. 1-2 

- 

Tm 

A. 1-3 

- 

SF.  =  SF.I 

1  Jlj=i  , 

SFj=2^"Sl+S2+S3"S4+S5+S6^ 

SFm  =  SFj 1 . 

J  lj=m 

S1 

■EDI 

- 

S2 

| 

warn 

S3 

A. 3-3 

warn 

S4 

mam 

warn 

S5 

A. 3-6 

A. 3-2 

S6 

HOI 

Ksa 

®i 

mi 

cd 

A. 1-4 

- 

EDx 

A. 1-5 

- 

®y 

A.  1-6 

- 

b. 

i 

i 

- 

mbm' 
— — 

c  =  speed  of  light  in  a  vacuum 
=  2.99792458xl08  m/sec 

n  =  surface  refractive  index 
=  1.000338 


ED 


=  13796.90  psec 
=  28094.50  psec 
=  41967.30  psec 


(A. 1-4) 
(A. 1-5) 
(A. 1-6) 
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Figure  A. 1-1  Mixed  Path  TD  Geometry 

JU 

SF^  =  Secondary  Phase  Delay  (SF)  from  1 
secondary  station- to-user  (nsec) 

SFm  =  Secondary  Phase  Delay  (SF)  from  master 
(m)  station- to-user  (psec) 

b.  =  TD  bias  associated  with  the  iC^  second- 
1  ary  station  (nsec) 


Note,  ED.  (i  =  w,  xory)  is  the  published  emission  delay  at 
the  i  secondary  station  (Ref.  1).  The  calibrated  values  of 
the  TD  biases,  b^ ,  for  the  GRB  and  LRB  models  are  given  in 
Table  A. 1-2. 


^Throughout  this  report,  SF  denotes  the  total  secondary  phase 
delay  of  the  groundwave  signal  propagating  over  any  land  and/or 
water  path. 


The  numerical  values  of  TD  biases  are  different  for 
each  model  and  each  TD  component  (i.e.,  TDW,  TDX  or  TDY).  In 
addition,  land  (i.e.,  all-land),  and  mixed  (i.e.,  land/sea  water), 
TD  paths  use  different  values  for  the  TD  bias,  as  tabulated  in 
Table  A. 1-2.  The  biases  presented  in  Table  A. 1-2  apply  when  both 
station  signal  paths  forming  a  TD  are  either  land  or  mixed. 


TABLE  A. 1-2 

CALIBRATED  TD  BIASES  FOR  GRB  AND  LRB  MODELS 


TD  BIAS 

-  jjsec 

TD 

rr\UDAXTP\TT 

STATION -TO- SITE  SIGNAL  PATH 

LAND 

MIXED 

GRB 

LRB 

GRB  . 

LRB 

TDW 

-0.854 

-0.394 

* 

* 

TDX 

-1.648 

-0.262 

-1.173 

-0.425 

TDY 

0.085 

0.294 

-0.353 

-0.470 

♦Available  to  estimate  mixed  path  TD  bias  for  TDW. 


Formulas  to  compute  a  mixed  path  SF  from  land  and  sea 
water  path  SFs  are  given  in  Section  A. 2.  The  calibrated  algo¬ 
rithms  for  the  SFs  of  the  GRB  and  LRB  models  are  presented  in 
Section  A. 3.  The  measured,  computed  and  residual  TDs  at  the 
TD  measurement  data  sites,  obtained  with  the  GRB  model,  are 
given  in  Section  A. 4. 


A. 2  MIXED  PATH  SF  EQUATION 

Using  Millington's  method  (Eqs.  2.2-1  through  2.2-3), 
the  formula  for  a  mixed  path  SF  from  the  (=  w,  x,  y  or  m) 
station  to  a  user  is 
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j  «  J  «i.<VV  *  SFL(TLJ,Pj)  -  SFL(Tsj) 


+  SFs(Tj>  +  SFS(TSJ>  -  SFs(Tlj)  psec  (A. 2-1) 


where  SF,(T,p)  is  the  SF  of  a  land  path  of  length  T  (Msec)  and 
bearing  angle  p  at  the  station;  and  SFg(T)  is  the  SF  of  a  sea 
water  path  of  length  T.  In  general,  the  functions  SF^  and  SFg 
may  be  station  specific.  The  path  lengths  Tt  . ,  Tc.,  T.  in  Eq. 
A. 2-1  refer  to  the  land  segment,  sea  segment  and  total  length, 
of  a  signal  path  from  the  station  to  a  user,  respectively. 

Section  A. 3  give*  formulas  to  compute  each  of  the  six  terms 
included  in  Eq.  A. 2-1. 


A. 3  EQUATIONS  OF  MIXED  PATH  SF  COMPONENTS 


A. 3.1  Term  S, 


The  term  S,  in  Eq.  A. 2-1  is  the  SF  of  a  land  path  of 

*  fU 

length  Tgj  (Msec)  from  the  j  station  and  is  given  by 


S1  =  SFL(TSj> 

=  Pf7-95  +  0.439  +  0.00245  Tc  .  Msec 


(A. 3-1) 


Note,  for  an  all- land  or  all-sea  water  path,  S^  is  zero. 


*Term  S^  in  Eq.  A. 2-1  does  not  depend  on  the  bearing  angle. 


A.  3.2  Terms  S2 >  S ^  and 

Terms  S2,  and  in  Eq.  A. 2-1  are  the  SFs  for  sea 
water  path  lengths  T . ,  Tgj  and  T^j ,  respectively.  The  sea  SF 
model  associated  with  both  the  GRB  and  LRB  models  depends  only 
on  range  to  the  station  in  the  following  manner: 


s2  =  SFS(T.) 

53  =  SFS<TSj> 

54  »  SFS<TLj> 


(A. 3-2) 
(A. 3-3) 
(A. 3-4) 


where  the  sea  water  SF  for  a  path  of  length  T  (psec)  is  given  by 


a_i 

-Tjr-  +  aQ  +  a^  T  psec,  if  10  <_  T  £  540  psec 


SFg(T)  = 


(A. 3-5) 


+  aQ  +  a^  T  psec,  if  T  >  540  psec 


The  calibrated  values  of  the  sea  coefficients  associated  with 
the  GRB  and  LRB  models  are  given  in  Table  A. 3-1.  Note,  for  an 
all -land  path,  the  terms  S2,  S^  and  S^  are  each  zero;  for  an 
all-sea  water  path,  however,  only  the  term  is  zero. 

A. 3. 3  Terms  andS^ 

Terms  and  in  Eq.  A. 2-1  are  the  SFs  of  land  paths 
of'  lengths  T.  and  Tt  . ,  respectively.  In  both  the  GRB  and  LRB 
models,  these  terms  include  bearing  angle  dependence  in  addi¬ 
tion  to  the  range  dependence.  Terms  S^  and  S^  are  given  below 
by  SFs  associated  with  a  land  path  of  -length  T^  (psec)  and 

subtending  a  bearing  angle  of  £  .  (measured  positive  clockwise 

th  J 

from  north)  at  the  j  station: 


TABLE  A. 3-1 

CALIBRATED  VALUES  OF  SEA  SF  MODEL  ‘ 
COEFFICIENTS  FOR  GRB  AND  LRB  MODELS 


COEFFICIENT 

MODEL 

GRB 

LRB 

a-l 

3.188 

2.885 

ao 

-0.594 

-0.387 

al 

0.000329 

0.000332 

f 

a-l 

128.8 

130.4 

1 

ao 

0.187 

-0.012 

1 

al 

0.000652 

0.000660 

55  =  SFL<Tjf  $.)  (A. 3-6) 

56  =  SFL(TLj’  Pj>  (A-3-7) 

where  the  formulas  for  the  SFt  of  a  land  path  of  length  T. 
(Msec)  and  bearing  p^  (deg)  for  the  GRB  and  LRB  models  are 
given  by 


GRB  Model: 


SFL(Tj,  Pj)  =  A0  +  A1  Tj  +  lcj£  sin  +  Dj£  cos  j  I  Msec 

(A. 3-8) 


LRB  Model: 


SF^T^pj)  =  Aq  +  [Ax  +  Bx  f.(Pj)J  Tj  psec  (A.  3-9) 


The  function  f.(p.)  is  zero  for  j  =  m,w  and  y.  However,  f . ( p  ) 
J  J  J  " 

is  unity  for  bearing  angles  between  3  and  150  deg  and  zero  for 


all  other  bearing  angles.  Calibrated  coefficient  values  of 
the  GRB  and  LRB  models,  Eqs.  A. 3-8  and  A. 3-9,  are  as  given  in 
Table  A. 3-2. 


TABLE  A. 3-2 

CALIBRATED  VALUES  OF  LAND  SF  MODEL 
COEFFICIENTS  FOR  GRB  AND  LRB  MODELS 


*Not  Applicable. 


tTheoretical  Values  (unobservable  from  the 
calibration  data). 


A. 4  COMPUTED  TDs  AT  CALIBRATION  DATA  SITES 

Tables  A. 4-1  and  A. 4-2  present  the  measured  and  com¬ 
puted  TDs,  and  TD  residuals  at  the  land  and  sea  calibration 
data  sites,  respectively.  The  computed  TDs  and  the  residuals 
were  obtained  with  the  GRB  model  which  has  been  selected  as 
the  West  Coast  Grid  calibration  model  (see  Chapter  4).  (Note 
Fig.  4.4-1  shows  plots  of  the  site  residuals  tabulated  in 
Tables  A. 4-1  and  A. 4-2.) 
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TABLE  A. 4-2 


MEASURED,  CALIBRATED  AND  RESIDUAL  TIME  DIFFERENCES 

AT  SEA  DATA  SITES 
(GRB  MODEL) 


T-3548 


TIME  DIFFERENCE  -  usee  | 

TASC 

SITE 

TDX 

TOY  | 

IDENTI¬ 

FICATION 

MEASURED 

CALIBRATED 

RESIDUAL 

MEASURED 

CALIBRATED 

RESIDUAL 

23 

27808*800 

27608.702 

0.096 

41916.000 

41915.664 

0.336 

2< 

27801.100 

27800.991 

0.109 

41920.000 

41919.625 

0.375 

25 

27913.400 

27913.421 

-0.021 

41698.300 

41696.198 

0.102 

25 

27835.700 

27635.474 

0.226 

41796.300 

41796.295 

0.005 

28 

27976.400 

27976.117 

0.263 

41592.400 

41592.217 

0.163 

29 

27 880.300 

27679  «  904 

0.396 

41699.000 

41698.699 

0.101 

31 

27965.400 

27965.311 

0.089 

41546.000 

41548.369 

-0.389 

32 

28036. 200 

28035.377 

0.823 

41410.500 

41410.958 

-0.458 

33 

27896.000 

27896.202 

-0.202 

41615.700 

41616.141 

-0.441 

34 

28094.900 

26094.947 

-0.047 

41294.500 

41294.593 

-0.093 

35 

28062.800 

28062.670 

0.130 

41339.100 

41339.212 

-0.112 

38 

28002.900 

28002.263 

0.637 

41419.400 

41419.384 

0.016 

37 

28165.400 

26 166.067 

-0.667 

41122.500 

41122.761 

-0.261 

38 

27976.600 

27976. 187 

0.413 

41450.900 

41450.946 

-0.046 

40 

28174.400 

26174.928 

-0.526 

410 36.000 

41036.355 

-0.355 

41 

26199.200 

26199.780 

-0.58Q 

40981.300 

40981 .494 

-0.194 

42 

26244. 100 

26244.715 

-0.615 

40846.100 

40845.359 

0.741 

44 

28273.500 

28273.992 

-0.492 

40675.300 

40674.968 

0.332 

45 

26120.400 

28120. 164 

0.236 

40975.600 

40975.866 

j 

-0.286 

48 

28271.900 

28272.205 

-0.305 

406 38.200 

40637.741 

0.459 

47 

28097.400 

26097.456 

-0.056 

40974.300 

40974.495 

-0.195 

48 

28262.800 

26262.762 

0.018 

40601.000 

40600.777 

0.223 

49 

28227.400 

28227.345 

0.055 

40674.000 

40674.043 

-0.043 
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